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Abstract. Binary mixtures of electrons with protons and of electrons with photons subjected to 
external electromagnetic fields are analyzed by using the Anderson and Witting model equation. 
The relativistic laws of Ohm and Fourier are determined as well as general expressions for the 
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relativistic completely degenerate electrons. 
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INTRODUCTION 

The analysis of non-relativistic and relativistic ionized gases by using the Boltzmann 
equation is a very difficult subject, since it refers to a system of coupled nonlinear 
integro-differential equations for the distribution functions. Simpler model equations 
for the collision term have been proposed in the literature in order to overcome the diffi- 
culties of the Boltzmann integro-differential equation. The model equations simplify the 
structure of the collision term but maintain its basic properties. For the non-relativistic 
Boltzmann equation the most widely known model is the BGK model which was for- 
mulated independently by Bhatnagar, Gross and Krook [1] and Welander [2]. The first 
extention of the non-relativistic BGK model to the relativistic case was proposed by 
Marie [3]. Although the non-relativistic limiting case of the Marie's model recovers the 
non-relativistic BGK model, in the case of particles with zero rest mass the relaxation 
time of the distribution function tends to infinity. This shortcoming was found by An- 
derson and Witting [4] who proposed a new model equation. 

In this work we follow [5] and analyze binary mixtures of electrons and protons 
and of electrons and photons subjected to external electromagnetic fields within the 
framework of Anderson and Witting model equation. These two systems are important in 
astrophysics since they could describe magnetic white dwarfs or cosmological fluids in 
the plasma period and in the radiation dominated period. By using the Chapman-Enskog 
methodology we determine Ohm and Fourier laws in the presence of electromagnetic 
fields and general expressions for the electrical and thermal conductivities for relativistic 
non-degenerated and degenerate binary mixtures of electrons with protons and electrons 
with photons. Furthermore, explicit expressions for these coefficients are given for 



the particular mixtures: (a) a non-relativistic mixture of protons and non-degenerate 
electrons; (b) an ultra-relativistic mixture of photons and non-degenerate electrons; (c) 
a non-relativistic mixture of protons and completely degenerate electrons; (d) an ultra- 
relativistic mixture of photons and completely degenerate electrons and (e) a mixture of 
non-relativistic protons and ultra-relativistic completely degenerate electrons. 



RELATIVISTIC UEHLING-UHLENBECK EQUATION 

Let us first consider a single relativistic quantum ideal gas in a Minkowski space charac- 
terized by metric tensor 7]"^ with signature diag(l,- 1,-1,-1). In the phase space spanned 
by the space-time coordinates (x") = (c?,jc") and momentum four- vector (p") = {p^,^) 
the state of the relativistic quantum gas is characterized by the one-particle distribution 
function f{x^ , p") = /(x, p, t) , since the length of the momentum four-vector is given by 
mc so that = -\/|pp+"m^?. The number of particles at time t in the volume element 
d^x about X and with momenta in the range d^p about p is given by /(x, p, t)d^xd^p. 

The space-time evolution of the one-particle distribution function /(x, p, = / in the 
phase space is given by the Boltzmann equation (see e.g. [6, 7]) 

^df dfK^ 



where m denotes the rest mass of the particle and is the Minkowski force which 
acts on the particles of the gas. Furthermore, Q is a term which takes into account 
the collisions of the particles. For a relativistic gas which obeys the classical statistical 
mechanics it is given by 

Q = I {f'j'-U)Fada^. (2) 
./ ^ p*o 

In the above equation we have introduced the abbreviations = /(x,p'^,r), /' = 
/(x,p',f), /* = /(XjPhc,?)' / = fi^^V:^): where p and p* denote the momenta of two 
particles before a binary collision and p' and p'^ are the corresponding momenta after col- 
lision. The pre and post coUisional momentum four-vectors are connected by the energy- 
momentum conservation law p" +p^ = p"^ + p'". Furthermore, F = \/ {p^PaY ~ fn^C^ 
is the invariant flux, which in the non-relativistic limiting case is proportional to the mod- 
ulus of the relative velocity. The differential cross-section and the element of solid angle 
that characterize the binary collision are denoted by o and dO., respectively. 

The collision term Q for a gas whose particles obey quantum statistics may be 
motivated as follows. First we note that the volume element in the phase space d^xd^p 
is a scalar invariant, but when quantum effects are taken into account in a semi-classical 
description, we divide the volume element by h^, where h = 6.626 x 10^-^^ J s is the 
Planck constant. Hence we write d^xd^p/h^, which is also a scalar invariant. The term 
d^xd^p I may be interpreted as the number of available states in the volume element 
d^xd^p. For particles with spin s there are more states, corresponding to the values that 
the spin component on a given axis can take and we have to introduce the degeneracy 



factor gs- Hence the number of available states is given by 



where 



g: 




2s +1 for 
2 s for 



m = 0. 



(3) 



In quantum mechanics a system of identical particles may be described by two kinds 
of particles: bosons and fermions. Bosons have integral spin, obey the Bose-Einstein 
statistics and include mesons (pion, kaon), photons, gluons and nuclei of even mass 
number like helium-4. Fermions have half-integral spin, obey the Fermi-Dirac statistics 
and include leptons (electron, muon, tau), baryons (neutron, proton) and nuclei of 
odd mass number like helium-3. The main difference between bosons and fermions in 
quantum statistical mechanics refers to the occupation number of a state. Any number 
of boson particles may occupy the same state, while fermion particles obey the Pauli 
exclusion principle and at most one particle may occupy each state. 

In order to incorporate the statistics of bosons and fermions into the collision term, we 
begin to analyze fermions and note that due the Pauli exclusion principle, the phase space 
is completely occupied if the number of the particles in d^xd^p is equal to the number of 
available states fd^xd^p = g^d^xd^p/h^, so that / = g^/h^ . Hence, (1 — fh? / g^) gives 
the number of vacant states in the phase space. If the number of particles that enter 
the volume element d^xd^p in phase space, as a consequence of a binary collision, is 
proportional to this quantity must be multiplied by the number of vacant states 
which is proportional to (1 — fh^ / g^){ \ — fjr'/gf). Hence the following substitution in 
the collision term of the Boltzmann equation must be consider: 



for the particles that leave the volume element d^xd^p in phase space. 

To include the apparent attraction between the boson particles - due to the statistics of 
indistinguishable particles with no restrictions on the occupation of a state - the factor 
(1 —fh^/gs) must be replaced by (l+//?^/^s) . Hence we can write from the Boltzmann 
equation (1) and from the above conclusions the relativistic Uehling-Uhlenbeck equation 
(for the non-relativistic Uehling-Uhlenbeck equation see [8]) 




(4) 



On the basis of the same reasoning we have to substitute 




(5) 




(6) 



where e is defined through 



+ 1 for Bose-Einstein statistics; 

— 1 for Fermi-Dirac statistics; 

for Maxwell-Boltzmann statistics. 



(7) 



At equilibrium the number of particles that enter and leave the volume element in the 
phase space must be equal to each other, so that the quantity within the brackets in (6) 

must vanish. Equivalently, In /(^^ /(I + £f^^^h^ /gs) must be a summational invariant 

- i.e., a function that obeys the relationship \]/ + \]/^ = \]/' + - where f^^^ denotes the 
equilibrium distribution function. For summational invariants there exists the following 
theorem (see e.g. [7]): A continuous and differentiable function of class W{p"') is a 
summational invariant if and only if it is given by ^^{p") — A + Bap'^, where A is an 
arbitrary scalar and 5« an arbitrary four- vector that do not depend on p". Hence we 
have 



In 



where a = —A — \n{gs/h^). 

For the determination of a and 5" we refer to [7]. Here we give only the results that 
a = ji/kT and 5" = U^/kT, where /i is the chemical potential, T the temperature, k the 
Boltzmann constant, and the four- velocity (with U"Ua = c^). Hence, the equilibrium 
distribution function reads 



/ 
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(9) 



when £ = and e = =i=l, respectively. The relativistic Maxwell-Boltzmann distribution 
function (9)i was obtained by Jiittner [9] in 1911 and the relativistic Fermi-Dirac (+) 
and Bose-Einstein (— ) distribution function (9)i was deduced by him [10] in 1928. 

The extension of the Uehling-Uhlenbeck equation to a mixture of r constituents is 
straightforward. We introduce an one-particle distribution function for each constituent 
of the mixture = /(x, Pa, {a — 1, . . . , r) which must satisfy the equation 
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(10) 



Above it was supposed that the external force that acts on the particles of electric charge 
q^ is of electromagnetic nature. In this case the Minkowski force reads 



rrin 



(11) 



where F"^ is the electromagnetic field tensor. 

Now we introduce the moments of the distribution function, which are the partial 

particle four- flow A/^ and the partial energy-momentum tensor Ta . They are defined 
through: 



J PaO J PaO 

The corresponding quantities for the mixture read 



(12) 



A^« = ^A^«, r«^ = ^r«/^. (13) 

a=l 0=1 

In the analysis of ionized gases it is also important to introduce the electric charge 
four- vector which is defined in terms of the partial particle four-flows A^" and of the 
partial electric charges as 

J^-t'laK- (14) 

fl=l 

The balance equations for the particle four-flow and of the energy-momentum tensor 
of the mixture are obtained by multiplying (10) by c and cp", respectively, and by 
summing the resulting equations, yielding 

^«iV« = 0, dpT^P = -F^P t = -F^^Jp- (15) 

^ a= 1 ^ 

Equation (15)i is the conservation law for the particle four-flow of the mixture. Equation 
(15)2 when compared with the balance equation for the energy-momentum tensor of the 

electromagnetic field has an opposite sign on its right-hand side. However, if we 
denote the energy-momentum tensor of (15)2 by an index pt - that refers to the particles 
- we get the conservation law (see Landau and Lifshitz [11]): 

da{Tpf+T,l^)=0, (16) 

which means that the sum of the energy-momentum tensors of the particles and of the 
electromagnetic field satisfies a conservation equation. 



LANDAU-LIFSHITZ DECOMPOSITION 

The decomposition of the partial particle four-flow and of the partial energy-momentum 
tensor proceeds by introducing the four- velocity U'^ and the projector A"^ defined by 

A^p ^ ^ap _ }_jja^l3 ^ such that A«^t/^=0. (17) 



In the Landau-Lifshitz description [12] the partial particle four-flow and the partial 
energy-momentum tensor may be decomposed according to 

< = n,f/«+7«-^, (18) 

Tf = Pi"^^ - {Pa + C^a)A-P + + hjS - ^^^) 

Above we have introduced the following quantities for the constituent a in the mixture: 

particle number density Ua, diffusion flux 7^, pressure deviator pi"^\ pressure pa, non- 
equilibrium pressure (Da, heat flux q", energy per particle ea and enthalpy per particle 
ha = ea + Pa/^a- The corresponding quantities for the mixture are given by the sums 



n=l^na, ;^<«^> = £pi"^\ P=l,Pa, QJ=£qJ., (20) 

a=l a=l a=\ a=l 

r r r 

ne=Y, naea, ila+haJa)i nh=Y, ^aha- (21) 

a=\ a=\ a=\ 

The sum of (18) and (19) over all constituents of the mixture lead to the following 
decompositions of the particle four-flow and energy-momentum tensor of the mixture 

A^« = „[/«-i!^, =p(«/5)_(p + Qj)A«/^-K^J[/«t//5, (22) 

h 

thanks to the constraint that there exist only {r — I) partial diffusion fluxes that are 
linearly independent for a mixture of r constituents, namely, 



I 7« = 0. (23) 



a=l 



We may also define the electric current four- vector in terms of the partial diffusion 
fluxes 7" and of the partial electric charges qa as 

/" = i;qA"- (24) 

We refer to the works of de Groot and Suttorp [13] and of van Erkelens and van 
Leeuwen [14] and decompose the electromagnetic field tensor F"^ into one part which 
is parallel to the four-velocity U"' and another which is perpendicular to it, i.e. 

= 1 [p'^yUyU^ -F^myU''^ +A«F^^A^. (25) 

Furthermore, by introducing the tensors and 5"^ defined by 

E'' = -F''PUii, = -A«F^^ , (26) 



we may write the electromagnetic field tensor as 

pap ^If^^ayP _^l5^a'^_gap ^27) 

If we consider a local Lorentz rest frame where ([/") = (c,0), equations (26) imply 
that 

(£«) = (0,E), B^"=B"^ = 0, B'j = -ce'j^Bk, (28) 

and we can identify E" with the electric field E and 5"^ with the magnetic flux induction 
B. 

Due to the fact that F"^ is an antisymmetric tensor F(^^U"'U^ = 0, it follows from 
(26) and (27) the relationships 

EaU^^O, B^pUl^^O, and = (29) 

CHAPMAN-ENSKOG METHOD 

Since we are interested to derive the laws of Fourier and Ohm for a binary mixture of 
electrons and protons and of electrons and photons, we have to made some simplifica- 
tions of our model, which are enumerated below: 

1 . the electric current four- vector (24) for a binary mixture of electrons {a = e) and 
protons {a = p) may be written as 

/" = -2e7f , (30) 

since the relationship between the diffusion fluxes reads = —J" and the elec- 
tric charges are given by = — e, q^ = e, with e denoting the elementary charge. 
Furthermore, we shall analyze the so-called Lorentzian plasma [15] where the col- 
lisions between the electrons may be neglected in comparison with the collisions 
between the electrons and protons. A Lorentzian plasma must fulfill the condition 
that the mass of one constituent is much larger than the mass of the other con- 
stituent. Here we have that mp/mg ^ 1836, where and denote the electron 
and proton masses, respectively. Moreover, we shall assume a locally neutral sys- 
tem where q^Ue + q^^p = 0, which implies that = rip-, 

2. the electric current four-vector (24) for a binary mixture of electrons (a = e) and 
photons (a = 7) , reduces to 

/" = -e7f, (31) 

due to the fact that the electric charge of the photons is zero (q^ = 0). Further- 
more, the collisions between electrons can also be neglected in comparison to the 
collisions between electrons and photons, which is the Compton scattering; 

3. the partial heat fluxes of the protons and of the photons are negligible in comparison 
with the partial heat flux of the electrons so that we can write from (21)2 that the 
heat flux of the mixture reduces to 



= ^" + (he - hh)J^, with b = p,r- 



(32) 



For simplicity we shall adopt the Anderson and Witting model equation [4] for the 
electrons instead of using the relativistic Uehling-Uhlenbeck equation (10). Hence, by 
taking into account the above considerations we write the space-time evolution of the 
distribution function for the electrons as 

adfe ^paH dfe _ U^'Pea . (0). 



where Xeb with b = p or b = y is the mean free time between collisions of electrons- 
protons or electrons-photons, respectively. In the above equation /j*^^ is the equilibrium 
distribution function of the electrons which reads 

/i'^ = ^^ ^ (34) 



by considering that the electrons obey the Fermi-Dirac statistics. Above, T denotes the 
temperature of the mixture, /i^ the chemical potential of the electrons and the factor 2 
refers to the degeneracy factor of the electrons. 

Once we know the equilibrium distribution function of the electrons we may calculate 
the values of the fields at equilibrium: particle number density Ug, energy density Ugeg 
and pressure pe defined by 



J PeO 



fpffi"''-^, (35) 



n,e, = ^UaUpTf = \uaUpc [ p^p^fP^. (36) 

C C J PeQ 

^ 3 J PeO 

The calculation proceeds as follows: we consider a local Lorentz rest system where 
= (c, 0) so that the particle number density of the electrons (35) reduces to 

ne = /^ exp(-^ + ^) + l 'P-''""^^^^'^^'P-'' ^''^ 

where we have introduced the spherical coordinates 0<\j/<n, 0<X^^^ and 
< I Pel < oo. Now we change the integration variable by introducing a new variable 
t> defined through 

Ipel = mgCsinhtJ, so that = ^^cosht?, (39) 

kT 

where = mgC^/kT is the ratio between the electron rest energy mgC^ and the thermal 
energy of the gas kT. When ^ 1 the electron behaves as a non-relativistic gas, while 
when ^ 1 it behaves as an ultra-relativistic gas. The change of variables and the 
integration of (38) in the angles x and i// leads to 

„ fmeC\^ sinh^tJcoshtJJt? Sn , ,3^ 



In the above equation we have introduced the electron chemical potential /i^ 
in units of kT and the integral Jnm{^e-, Me ) defined by 



°° sinh"-i^cosh'"-[?Jt? 



/o exp(-jU* + C,e cosh i^) + 1 
Following the same methodology we get that 



(41) 



8;r 



(42) 



Now we shall determine from (33) the non-equilibrium distribution function for the 
electrons by adopting the Chapman-Enskog methodology. For that purpose we search 
for a solution of the form 

(43) 



where the deviation from the equilibrium distribution function is considered to be a small 
quantity, i.e., I^^l ^ 1- If we insert (43) into the Boltzmann equation (33) we get 



(0) 



Pep 
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(0) 



(44) 



where we have not taken into account the term d^ejdx'^, since it is not our aim in 
deriving constitutive equations which are functions of second-order derivatives (Burnett 
equations). The above equation can be written as 
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<^e. (45) 



where we have not considered the term E'^d^g/ dp'^ , since it refers also to a second- 
order term. Furthermore, we have introduced in the above equation the electron cy- 
clotron frequency (Oe = eB/nie - where B is the modulus of the magnetic flux induction 
- and the differential operators D = U"da and V" = A^^dp . 

In this work we are interested in the derivation of the laws of Fourier and Ohm, so 
that we can restrict ourselves to the thermodynamic forces that are four- vectors, namely 



v«r: 



and 



kT 

£«__V« 
e 



(El. 
\kT 



(46) 



the first being a combination of a temperature gradient and an acceleration, while 
the second refers to a combination of an external electric field and a gradient of the 



chemical potential of the electrons. Hence, we obtain from (45) that the deviation from 
the distribution function may be written as 



(47) 



Up to terms in (cOeTeb/B)^ the four-vector A" is given by 
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mgC / COeTeb 
U^PeS^ B 



(48) 



Equation (47) together with (48) represent the deviation of the distribution function 
of the electrons as a function of thermodynamic forces that are four-vectors. We shall 
use the distribution function (43) in the following section in order to determine the laws 
of Ohm and Fourier. 



OHM AND FOURIER LAWS 

The determination of the diffusion flux 7" and the heat flux of the electrons proceeds 
by noting that (12), (18) and (19) lead to 




The insertion of the distribution function of the electrons (43) together with (47) and 
(48) into (49) and (50) and integration of the resulting equations, implies a system of 
equations for 7" and which is used to determine the heat flux of the mixture (32) 
and the electric current four- vector (30) or (31). From this system of equations it follows 
Fourier and Ohm laws 

= A«/^v^r+T«%^, /« = a«%^ +a«^v^r, (5i) 

respectively. Above A"^ is a tensor associated with the thermal conductivity, a"^ is 
the electrical conductivity tensor, while the tensors T"^ and ^2"^ are related with cross 
effects. We may represent the general expressions for the above mentioned tensors as 



+ {a3,Z73,C3,J3}5"^5?, (52) 



where the scalar coefficients ai through are given below: 
1 . Coefficients associated with A"^ 



^ 3h^hbkT^ ■ ^^^^ 



2. Coefficients associated with T"^ 



SnmyTebho { . hb , 



3h^hbkT \ MgC^ 
STtm^c^Tebhe f , hb . \ /(OeT^z,\ 

= - 3/^3 w (57) 



3h^hbkT y^-^ " ) ■ ^^^^ 



3. Coefficients associated with cr"^ 

%nrnlc^Xebne{Z+\)Q^ ( fibhb . \ 

- 3h^nhbkT y^'^^'^w^^'-')' ^^^^ 

_ Snm'lc'^T;^bne{Z+\)e^ f , ribhb , \ ( 0)eT:eb\ 

'2 - 3^wr 1,^4-1 +;^:;;^:?^4-2j(^^j, (60) 

^ 8;rm^c^Tej,ne(Z+ l)e^ / . ^ n^hb ^, \ /W^y ^^^^ 
3h^nhbkT \ netngC^ ^ cB J 

4. Coefficients associated with ^2"^ 

^1 ^ l^^4i + ;^:/4oJ, (62) 

- 3)^Wr^ (^^40 + ;^J4-iJ(^^j, (63) 

8;r7«^c'^Tefone(Z+l)e / n^/z^ \ /COeT^^X^ 

^3 = M^W^ 1,^4-1 + ;^:/4-2j(^^j • (64) 

In the above equations J7J*„ represents the partial derivative of (41) with respect to the 
chemical potential of the electrons /i* = jig/ (kT) in units of kT. Furthermore, we have 
introduced the abbreviation Q = mgC^ /{kT) which refers to the ratio between the rest 
energy of the electrons m^c^ and the thermal energy of the mixture kT . We note that 
in all above equations one has to consider Z = 1 for binary mixtures of electrons and 
protons and Z = for binary mixtures of electrons and photons. 



The thermal conductivity tensor X*^^ is obtained by eliminating from (51)i 
through the use of (51)2 by assuming that there is no electric current. Hence, we get 
a relationship between S'^ and V"T from (5 1)2 which may be used to write Fourier law 
as 

q" = X''I^VpT, where = eiT]«^ +^25"^ +e35«^5?. (65) 

Up to terms in [cOeteb / (cB)]^ the scalar coefficients ei through read 

aici—bidi a2c\ — bi{cid2 — C2di)—b2Cid\ ^^^^ 
e\ = , 62 = 5 , (66) 

a3c] - bi [di {c\ - C1C3) - ciC2d2] - c\{bidj, + bj,dx) - Cib2{c\d2 - C2d\) 
e-i = 5 • (67) 

In order to get a better physical interpretation of the components of the tensors, it is 
usual in the theory of ionized gases to decompose the thermodynamic forces V^T and 
Ea into parts parallel, perpendicular and transverse to the magnetic flux induction. To 
achieve this goal we follow van Erkelens and van Leeuwen [14] and introduce the dual 
B^^ of the magnetic flux induction tensor B'^^ defined by 

5«/3 _ l£«/5r55^^ (68) 

One may easy verify from (68) and (28) that in a local Lorentz rest frame the only 
non-zero components of 5"^ are 5^' = c5' since 5'' = and B^'^ = 0. 

The desired decomposition of the thermodynamic forces into parallel V"T, S^; 
perpendicular V" T, and transverse Vf T, parts read 

^ir - jr^f''^^''^ '--j^f'^^^'^ '''' 

= ^ r5"^^/3, (70) 

where J'" is an abbreviation for £" or V^T. In a local Lorentz rest frame (69) and (70) 
reduce to 

^0 = j-o = j-o^O, ^|i = -l(B-^)B, (71) 
^± = ;^[(B ■ ^)B - (B ■ B)^], = X B), (72) 



thanks to the relationship ^ B'y^By§/l = ca/B B = cB. From the above equations it is 

easy to verify that J^y is parallel to the magnetic flux induction B, J^x perpendicular to 
it while J-'t is perpendicular to both J^y and J^^. 
Now by using the following relationship 

(c5)27]«^ = (73) 



the Fourier and Ohm laws can be rewritten in terms of ^-"11^, J-"" and J-"". In fact, if we 

substitute (73) into the Ohm's law (51)2 and Fourier's law (65) 1 and make use of the 
definitions (69) and (70), it follows that the electric current four- vector and the heat flux 
can be written, without the cross-effects terms, as 

/« = a||£« + a^£« + a,£«, q"" = X^W^T + XiWlT + Xt^fT , (74) 

respectively. In the above equations the scalars are called the parallel, perpendicular and 
transverse components of the tensors, and their expressions are given by 

ra||=ci, aj_ =ci -C3(c5)2, Ot = C2{cB), 

\A||=ei, X^ = ei-e^{cBf, Xt = e2{cB). ^ ^ 

From the above formulas we shall obtain the parallel, perpendicular and transverse 
components of the electrical and thermal conductivities for binary mixtures of electrons 
and protons and of electrons and photons. 



ELECTRICAL AND THERMAL CONDUCTIVITIES 
Non-degenerate electrons 

Here we shall analyze two important cases, namely: a non-relativistic mixture of 
protons and non-degenerate electrons and an ultra-relativistic mixture of photons and 
non-degenerate electrons. We note that the chemical potential of the electrons in the 
non-degenerate case must fulfill the condition that e^'^^ ^ 1. 

1. A non-relativistic mixture of electrons and protons is identified by two conditions 
mp/nie ^ 1 and = nieC^ / (kT) ^ 1. In this case the transport coefficients read 
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(76) 
(77) 
(78) 
(79) 



The first relativistic corrections to the transport coefficients are related to the term 
1/ C,e and if fix our attention to the leading terms without the relativistic corrections, 
the electrical conductivities can be written from (76) and (77) as: 



2 

^ ^ep^e 



Ot = O\i{(0eTep), ^.L 



l + {COeTep)^'' 



(80) 



since we have considered COetep <^ 1. The expressions for the electrical conduc- 
tivities (80) are well-known in the theory of non-degenerate and non-relativistic 
ionized gases (see, for example, Cap [16]) and show their dependence on the mag- 
netic flux induction B through the electron cyclotron frequency £0^. Furthermore, 
the thermal conductivities (78) and (79) without the relativistic corrections become 

Note that the expression for the parallel thermal conductivity is well-known in the 
theory of non-relativistic gases which follow from a BGK model equation. 
2. An ultra-relativistic mixture of photons and non-degenerate electrons is character- 
ized by the condition = m^c^/ {kT) <^ 1. Here the transport coefficients reduce 
to 

^11=^^ = UnkT ' Uni^ ' ^^^^ 

h-^^- 3n, + 4n/ " (3n, + 4n,)2 " ^^^^ 

We infer from these equations that the parallel and perpendicular electrical and 
thermal conductivities coincide, while the transverse electrical and thermal con- 
ductivities are small quantities since they are proportional to ^e- 



Completely degenerate electrons 



All thermal conductivities vanish in the limit of completely degenerate electrons, since 
this behavior is connected with the well-known result from statistical mechanics that the 
heat capacity of a completely degenerate gas vanishes. For the electrical conductivities 
there exist three important cases to be analyzed which are: a non-relativistic mixture 
of protons and completely degenerate electrons; an ultra-relativistic mixture of pho- 
tons and completely degenerate electrons and a mixture of non-relativistic protons and 
ultra-relativistic completely degenerate electrons. We proceed to analyze the electrical 
conductivities for these cases. 

1 . A non-relativistic mixture of protons and completely degenerate electrons is iden- 
tified by ^ 1 and pf <^ mgC, where pf denotes the Fermi momentum of the 
electrons. Here we have 



Zne^Xeppl pI \ 



^ _ STte^TeppljcOeTep 



_ %ne^Xepp\ 



1 - 



(Oe^ep) 




Let us fix our attention to the leading terms of the electrical conductivities 



"^1 = ^;^^^' c., = a||«».T,„), a^^^-^—y, (86) 

since the term pf/{meC^) is a small quantity and the condition (Oe'T^ep ^ 1 holds. 
These equations show the dependence of the electrical conductivities on the mag- 
netic flux induction B through the electron cyclotron frequency (Og. 
An ultra-relativistic mixture of photons and completely degenerate electrons is 
characterized by the conditions <^ 1 and pf ^ rngC, and the electrical conduc- 
tivities for this case read 

" llnkTh^ \ UgCpF ) ' 

%Tle^Xeyne^eC{G)eTey)PF ( ^ '^kTny \ ^^^^ 



\2nh^ \ UgCpf 

We infer from the above equations that the parallel and perpendicular electrical 
conductivities are equal to each other, while the transverse electrical conductivity 
is a small quantity since it is proportional to ^g. 

A mixture of non-relativistic protons and ultra-relativistic completely degenerate 
electrons is also an important case since it could describe a white dwarf star. Here 
the conditions nip /me ^ 1 and pf ^ mgC hold and the electrical conductivities 
become 

showing that the parallel and perpendicular conductivities coincide and that the 
transverse conductivity is a small quantity, since it is proportional to ^ 1 . 



APPENDIX: INTEGRALS J' 



nm 



1. Non-degenerate case 

In this case e^f^" ^ 1 so that the integrals J^n reduce to 



rOO 

JnmiCe.^e) = / siuh" COSh'" I^cIt}, 

Jo 



(90) 



and the integrals J7j*„(^e,/i*) = J*^ can be expressed in terms of the modified 
Bessel functions of the second kind K„ ( Ce) = Kn and of their integrals Ki„ ( ^g) =Ki„ 
(see Abramowitz and Stegun [17] pages 376 and 483) as follows: 

J'l = 2^{K4- J*Q = -^{K^-Ki), J4_i = -^(i^2 -Ki2), (91) 
J:.2 = ^ (i^i + Kii - Ki3 ) , J4*_3 = ^ (Ki2 - Ki4) . (92) 

be be 



Moreover, the chemical potential of the electrons is given by 



%nmlckTK2 



(93) 



2. Completely degenerate case 

The integrals Jnm for this case reduce to 







sinh" cosh'" ^d^, with = arcosh W 1 + 



PF 
lUpC 



(94) 



where /'f is the Fermi momentum of the electrons. The integrals J7J*„ read 
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